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Abstract 



The pion-deuteron scattering length is computed to next-to-next-to-leading or- 
der in baryon chiral perturbation theory. A modified power-counting is then formu- 
lated which properly accounts for infrared enhancements engendered by the large 
size of the deuteron, as compared to the pion Compton wavelength. We use the 
precise experimental value of the real part of the pion-deuteron scattering length 
determined from the decay of pionic deuterium, together with constraints on pion- 
nucleon scattering lengths from the decay of pionic hydrogen, to extract the isovector 
and isoscalar S-wave pion-nucleon scattering lengths, a~ and a + , respectively. We 
find a' = [0.0918(±0.0013) ]M~ l and a + = [-0.0034(±0.0007) ]M~ l . 



PACS nos.: 13.75. Gx , 12.39.Fe 



1 Introduction 



The pion-nucleon (-7T-N) scattering lengths are quantities of fundamental importance in 
hadronic physics. They provide an important test of QCD symmetries and the pattern 
of their breaking. They also provide a crucial constraint on the 7r-N interaction and, 
as such, have an impact on our understanding of nucleon-nucleon (NN) scattering, the 
three-nucleon force, and pion-nucleus scattering. 

In the limit of exact isospin symmetry the threshold 7r-N amplitude may be written as 



7%, = 4tt(1 + //) [5 ba a + + ie bac r c a-\ , (1) 

where /i = M n /m ~ 1/7 is the genuine small threshold parameter, and a + and a~ are 
the isoscalar and isovector S-wave scattering lengths, respectively. While these quantities 
cannot be measured directly in 7r-N scattering experiments, the well-known effective field 
theory (EFT), baryon chiral perturbation theory (%PT), offers an ideal tool for extrapo- 
lating the 7r-N amplitude to threshold. Working up to fourth order in the dual momentum 
and pion-mass expansion, the values shown in Table [l] have been obtained from various 
analyses of 7T-N data |Q |2|, |3|, |[]. (Note that there are no fourth-order corrections to a~ ||. 
The difference between the 0(p 3 ) and 0(p A ) values in Table [I] arises from a refitting of 
low-energy constants at lower orders, as discussed in Ref. ||.) 





0(p) 


0(p 2 ) 


0(p 3 ) 


0( P 4 ) 


Fit 1 


0.0 


0.0046 


-0.0100 


-0.0096 


a + Fit 2 


0.0 


0.0024 


0.0049 


0.0045 


Fit 3 


0.0 


0.0101 


0.0014 


0.0027 


Fit 1 


0.0790 


0.0790 


0.0905 


0.0903 


a- Fit 2 


0.0790 


0.0790 


0.0772 


0.0771 


Fit 3 


0.0790 


0.0790 


0.0870 


0.0867 



Table 1: Convergence of the S-wave scattering lengths for x?T fits to the Karlsruhe ||, 
Matsinos 0, and VPI || phase-shift analyses of 7T-N data. 0(p n ) means that all terms 
up-to-and-including order n were included in the xPT fit. Units are M" 1 . Table taken 
from Ref. ||. 



An independent approach to the 7T-N scattering lengths involves analyzing pionic-atom 
level shifts and widths. In the Coulombic 7r~-p system, the strong-interaction shift in 
the energy can be used to infer a value for or. The recent Neuchatel-PSI-ETHZ (NPE) 
experiment finds: 

a - = [0.0905(±0.0042)]M^ 1 . (2) 



1 



Note that the error here is significantly smaller than the spread of values obtained from 
the analyses of 7T-N scattering (see Table [I]). Only the lifetime (or, equivalently, the width) 
of pionic hydrogen is sensitive to a + . The NPE experiment M finds the value: 

a+ = [-0.0022(±0.0043)]M 7 7 1 . (3) 

Thus, neither 7T-N scattering nor the 7r~-p atom provide a strong constraint on a + . The 
isoscalar scattering length may well be probed more directly in the 7r~-deuterium atom. 
The NPE measurement [n| of the pionic-deuterium atomic-level shift yields: 



a^ d = [-0.0261(±0.0005) + i 0.0063(±0.0007) JM" 1 (4) 

for the pion-deuteron (7r-d) scattering length — a measurement that is remarkably accurate 
given the usual level of precision in hadronic-physics experiments []• In order to precisely 
relate this number to a + an EFT of threshold n-d scattering is required. 

There has been remarkable recent progress in developing the EFTs relevant to multi- 



nucleon systems [|T|, [13|]. Among the advantages of this formalism is a quantitative method 
for estimating theoretical errors and a unified field-theoretic treatment of processes involv- 
ing different numbers of nucleons. The non-perturbative effects responsible for deuteron 
binding are accounted for in the EFT power-counting, with a meaningful quantification 
of the theoretical error. The multi-nucleon EFT relevant to momentum transfers of order 
the pion mass is in most straightforward generalization of single-nucleon baryon 

%PT |TJ], [T5|, |HJ. NN phase shifts and deuteron properties have now been computed 



to next-to- next-to-leading order (NNLO=0(p 3 ) for the NN potential, where p denotes a 
small momentum/pion mass) in this expansion O, [L^, |I9 |. Processes with external pions 
and photons have also been extensively studied [IB|. 

Consider 7r-d scattering at threshold [^, pT| . Two-pion-four-nucleon operators which 
contribute to 7r-d scattering enter at high orders in the EFT expansion and are therefore 
highly suppressed |2IJ . That these contributions are very small is crucial to the predictive 



power of the EFT, as the coefficients which determine the strength of these operators 
are unknown, and need therefore be determined from a separate nuclear observable, or, 
at some time in the distant futurity, from lattice QCD ||16||. This suppression of local 



operators has an associated benefit: the 7r-d scattering length is sensitive to low-energy 
constants which contribute to the S-wave 7r-N scattering lengths, a~ and a + . 

While baryon provides a rigorous theory of low-energy 7r-d scattering, there is 



an important kinematical subtlety in the threshold region; this leads to a new EFT f22 | 
which distinguishes between the low scales M v and 7, where 7 is the deuteron binding 
momentum. This straightforward generalization of baryon xPT turns out to be an ex- 
tremely useful tool in developing a precise theoretical relation between the 7r-d and 7r-N 
S-wave scattering lengths. 



1 For a general introduction to the 7r-d scattering length and its significance see Ref. |Tl| 
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Some additional remarks are in order. First, we should mention that the EFTs used 
in this paper can be applied to scattering and bound-state observables. For the direct 
extraction of the S-wave scattering lengths from pionic atoms one can also use the bound- 
state EFT as applied, for example, to 7r - -p atoms in Ref. [p3fl . Since that approach is 
designed to deal with the pionic-hydrogen system, it is expected to give better control of 
the theoretical errors. However, for the case of pionic deuterium one faces a genuine three- 
body problem, which so far has not been formulated in the bound-state EFT, whereas it 
is straightforward to deal with in our approach. Second, throughout this work, we neglect 
isospin violation. Isospin violation has been worked out in the 7r-N system to 0(p 3 ) |23J 



and in NN scattering to NLO p5|, [26], p7| but it goes beyond the scope of the present 



paper to systematically include it in the three-body system. 

In Section [2] we extend the baryon x?T analysis of the 7r-d scattering length to next-to- 
next-to leading order {0{p A )). We then formulate a modified power-counting in Section^ 
with which we derive an expression to extract the S-wave 7r-N scattering lengths from 
pionic-atom data. A comparison with other extractions of these quantities from the 
pionic-atom data is presented in Section |J We conclude in Section j|. 

2 The 7r-d Scattering Length in Baryon 



We decompose the 7r-d scattering length as pOL 



a-TTd — t^- ^yr>) ^ a " !Tn anp ^ ~^ a ( 6oos ') ^ a (3-feodj,) + i Im dird ■ (5) 

Each of the various contributions on the right-hand side of Eq. (|5|) will be discussed in 
turn. 

2.1 The isoscalar scattering length contribution 

We first consider the "two-body" contributions where the pion interacts with a single nu- 
cleon, while the other nucleon acts as a spectator. This contribution to a n d is proportional 
to the isoscalar S-wave 7r-N scattering length, 

a-Kn + a-Kp = 2a + . (6) 



This quantity has been computed to 0(p ) |28] and 0(p ) in baryon xPT. Charac- 



teristic diagrams are shown in Fig. [l|. The next-to- leading order (NLO= 0(p 3 ) for the 



amplitude) result is |28 
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where A = — 4ci + 2c2 + 2C3 and the Cj are low-energy constants. At this order there is a 
rather large spread in the values of a + resulting from fitting the low-energy constants to 
pion-nucleon phase-shift analyses |1|, 0, |3|, |J (see Table [I]). As pointed out above, mea- 
surement of the width of pionic hydrogen does not ameliorate this problem significantly, 
since the resulting value of a + is consistent with zero, and the error bar is too large to 
facilitate an accurate extraction of A. 

The pionic-hydrogen measurement does constrain a~ quite strongly. On the theoretical 



front, we have the 0(p ) expression for the isovector S-wave 7r-N scattering length || |29] 



M w AMI ( ^ g\\ Ml 



2F2 Fl \ 32m 2 / 16vr 2 F4 v ' 

TV TV \ / TV 

where D = d\ + d>2 + d% + 2d^ and the di are low-energy constants that have been fit to 
7r-N phase shifts in Refs. Jl| and ||. The NPE pionic-hydrogen measurement gives the 
value quoted in Eq. (Q), which allows for a determination of the combination of constants 
D. 



I H I 

(a) (b) 

\ 

I II I 

(c) (d) 

Figure 1: Characteristic two-body Feynman graphs which 
contribute to the 7i-d scattering length at 0{p 2 ) (a), 0(p 3 ) (b) 
and 0(p A ) (c),(d). The dots are leading-order vertices from 
£^jv, the shaded blobs are from C^n\ an d the solid square is 



from C^jr. 



2.2 The boost correction 

We now turn to the other contributions to a^d that can be reliably computed in baryon 
%PT. The piece u a^ oost ^" (sometimes known as the "Fermi motion" correction) is a con- 
tribution to a n d from the graph Fig. |](a), where the sliced blob represents a momentum- 
dependent vertex from the 0(p 2 ) chiral Lagrangian, C^jj. This vertex would vanish if the 
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7r-N system were at rest, but it contributes to 7r-d scattering, since the 7r-N amplitude 
must be boosted from the 7r-N center-of-mass frame to the 7r-d center-of-mass frame. This 
boost produces an 0(p A ) effect: 

where q = p—p' and p and p' are, respectively, the initial and final-state relative momenta 
of the two nucleons, and 

(d(q)U = J d 3 pd 3 p' ^ d (p) d{q) * d (p') , (10) 

where ^ is the deuteron wavefunction. The three-dimensional delta-function in Eq. (|9|) 
appears inside the expectation value since this boost effect is a "two-body" contribution 
to a nd . Note that here the momentum-space deuteron wavef unctions are normalized so 
that: 




2.3 Three-body effects 



I M M l 
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(b) 



(c) 



Figure 2: Three-body Feynman graphs which contribute to 
the 7r-d scattering length at 0(p 3 ) in baryon xPT. 

The dominant piece of the "three-body" amplitude, which occurs at 0(p 3 ) and is shown 
in Fig. |2|, was written down some time ago by Weinberg [^0 , 
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Nominal 0(p 4 ) corrections to these three-body effects are shown in Fig 
vertices appearing in these graphs are from C^ N 
according to the following pattern: 
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Amusingly, these amplitudes cancel 
= 0; (14) 
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Figure 3: Three-body Feynman graphs which contribute to 
the 7r-d scattering length at 0{p A ) in baryon x?T. The double 
lines indicate graphs evaluated in time-ordered perturbation 
theory. 
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M {a) +M {b) = 0; (15) 

M (f) + M ( 9 ) +M (k) +M(fi = 0; (16) 

M( h) +Mw+M® + M( k m) = 0\ (17) 

M H + M {o} = 0. (18) 

Thus there is no three-body correction to a na > at order 0(p A ). The sum of all corrections 
vanishes for a variety of reasons, among them the threshold kinematics and the isoscalar 
character of the deuteron. 



2.4 Results 

This completes the computation of a nc [ to 0(p 4 ) in baryon %PT Q The pieces not propor- 
tional to a + are a^body), which is simply the 0(p 3 ) expression given in Eqs. dl2|)-(|l3D, 
and a^oost), which is obtained by taking an expectation value of the one-body operator 
p 2 between deuteron wavefunctions — see Eq. @. The deuteron wavefunction computed 
with the NLO (0(p 2 )) %PT NN potential gives a consistent determination of a n d to 0(p 4 ). 
We have also shown results for the "NNLO*" x?T NN potential which yields a node-less 
deuteron wavefunction, and the NNLO NN potential, whose deuteron wavefunction 
has two nodes. These nodes signal the presence of unphysical, deeply-bound states. The 
numerical results for the various contributions to a n d are shown in Table |2|. The parame- 
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Table 2: Numerical results for the contributions to a^d given in Eqs. fl9|), (|12|) and flTB]) as 
well as other overlaps discussed in the text. All scattering lengths are in units of M~ x . A 
is the cutoff in the NN system. 



2 We relegate discussion of the imaginary part of a nc i to Section ^. 
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ter set = 92.4 MeV, g& = 1.2843 (note that we adjust gA from the Goldberger-Treiman 
relation), M w = M w + = 139.57 MeV and m = (m n + m p )/2 = 938.92 MeV is used in all 
expressions. For each of the three different ;\TT NN potentials we vary A, the cutoff in 
the NN system, in order to estimate our theoretical error. The three values (Fits 1-3) 
quoted for a(b oos t) correspond to values of c 2 fit to various phase shift analyses of 7T-N 
scattering || (Fits 1-3 in Table [I]). For comparison, in Table |3] we give numerical results 
for the same quantities computed using various "realistic" potential models. 
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Table 3: Numerical results for the contributions to a^d given in Eqs. (|9|),([T2]) and ( |i~3|) as 
well as other overlaps discussed in the text, for various potential models: Nijmegen |30| ; 
CD-Bonn 2000 M\; AV-18 [62. All scattering lengths are in units of Mr 1 . 



The following points are worth noting: 

1. The nodes in the NNLO wavefunction cause the anomalous change in a^ oost ) (by an 
order of magnitude) between the NLO and NNLO evaluations. This anomaly can be 
understood as follows. The boost correction to a na > is proportional to (p 2 ) w {. While 
the sum of (p 2 ) w f/m and (V) w f — the expectation value of the NN potential — is con- 
strained to be 2.225 MeV, the NNLO potential is very deep because it supports 
unphysical NN states bound by over one GeV. It naively follows that (p 2 ) w f will be 
much larger for the NNLO wavefunction than for wavefunctions without unphysical, 
deeply-bound NN states. Of course if renormalization is done properly, this sensitiv- 
ity to short- distance physics should be absent in observables up to subleading effects 
in the EFT. Given this technical issue, which is currently under investigation, in 
this paper we discard results obtained with the NNLO wavefunction. 

2. When evaluated with the chiral NLO and NNLO* wavefunctions the smaller pieces 
of the 0(p 3 ) three-body contribution, a 2 ^1 hodj ^, are rather different from the results 
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given in Table 1 of Ref. plj . In contrast, the difference between the evaluations 
of Table [| and those of Ref. is small and arises from the different value of 
qa used in that paper. The spread in a(|i body ) values in Tables || and reflects 
the sensitivity of this contribution to the (unphysical) short-range behavior of the 
deuteron wavefunction. 

Lastly, and most importantly, in all cases o^_ body ^ is much larger than any other 
contribution to — including the contribution from a + . This suggests that the 
baryon scaling of operators is not properly accounting for the characteristic 

scales present in the deuteron wavefunction. We will return to this important point 
in the next section. 



3 A Modified Power-Counting 

As pointed out by Weinberg , Feynman graphs of the type shown in Fig. 0(a) are larger 



than their naive baryon x?T scaling suggests because of infrared enhancements resulting 
from the anomalously small deuteron binding energy. For all of the NN wavefunctions 
considered in Table 0, a 2 ^_ body ^ is larger than a?|f. 6otJ \ by two orders of magnitude. This 
is mainly because in low-energy 7r-d scattering pions carry energy and momentum that 
puts external pion legs, as well as internal pion legs in diagrams such as Fig. 0(a), on or 
near the pion mass-shell. The amplitude in Fig. 0(a) would be infrared divergent if the 
deuteron were a zero-energy bound state. As it is, this divergence is softened by the fact 
that the two nucleons in the deuteron are off their mass-shell. But the amount by which 
they are off their mass-shell is far less than the pion mass, and therefore if these effects are 
to be correctly accounted for in the EFT, we must find a way to incorporate the infrared 
enhancement into the theory. 

The completely consistent way to do this is to construct an EFT where non-relativistic 
nucleons interact via contact operators, and pions near their mass-shell are included as 
massless excitations. Thus, in 7r-d scattering the nucleons are largely static and the 
successive scattering of a pion off the two nucleons generates a Coulomb-like interaction 



between them ||22|| . In order to have the incoming and outgoing nucleons as well as all 
the pions in a process near their respective mass-shells, there must be an equal number of 
pions entering and leaving each interaction vertex. Processes with interaction vertices that 
have a different number of pions entering than leaving are represented by local operators 
(on the scale of M^) in this "EFT with heavy pions" (FT/rEFT). H7rEFT is organized 
as an expansion in / ~f/M. K , since physics at the scale of the pion mass is "integrated 
out" and only appears implicitly in the theory via constants in the Lagrangian. Effects 
contributing explicitly to 7r-d scattering in this EFT are two-body 7r-N scattering and 



'three-body" effects in the multiple-scattering series for the nNN system [22]. Other 
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effects, such as ct(f l^y) ; are non-dynamical and are represented by a local two-pion- 
four-nucleon operator. This operator appears at leading order in the expression for a nd 
and its coefficient is undetermined in HzrEFT, since it is set by pion-scale physics. This 
counterterm fixes the scale dependence of a logarithm that is formally divergent in the 
infrared limit, 7^0. 

Here we will employ an economical way of recovering the hierarchy of scales 7 <C M n 
that motivates the construction of the H7rEFT. We will work in the pionful theory with 
a modified power counting that accounts for the infrared enhancement of graphs in the 
multiple-scattering series for the irNN system. In practice this means that we are able 
to fix the strength of the counterterm by matching the results of HzrEFT developed in 
Ref. [^] to those of baryon derived above. 

If we denote a generic momentum inside deuterium by q, then the baryon x?T power- 
counting is predicated upon | q\ ~ M n . One might instead expect | q\ ~ ^—B d m = 7, 
where B d and 7 are the deuteron binding energy and momentum, respectively. These 
two scalings would be compatible if the deuteron binding energy were of natural size. 
However, B d = —2.225 MeV is unnaturally small on the scale of hadronic physics, as 
is 7 = 45.7025 MeV. The true scaling is thus | q\ ~ M n /3, and this additional 1/3 
suppression has a dramatic effect on the ordering of operators. 



• MM* 

(a) \ (b) \ (c) \ 

Figure 4: Infrared-enhanced Feynman graphs which con- 
tribute to the 7r-d scattering length at 0(p 5 ) in baryon 
and 0(Q 3 ) in the modified power-counting. The small shaded 
blobs are vertices taken from £^y and the large shaded blobs 
are vertices taken from together with 0(p 3 ) one-loop 
graphs. 

Hence we will adopt a modified scaling of operators in which \ q\ ~ p 2 where p is the 
baryon xPT expansion parameter []. Amplitudes involving M n in the denominator are 
now suppressed, just as is the case in H7rEFT. The expansion parameter for this modified 
power-counting is denoted by Q, where Q is the same size as p in baryon but, of 
course, leads to a different ordering of operators. 

3 Yet another power-counting for 7r-d scattering in which pions are purely perturbative has been studied 
in Ref. [p3[ . Unfortunately, it is now known that pions are strongly non-perturbative in the deuteron 
channel |34|. 
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It is easy to check that in the modified power-counting a?^_ body ^ is 0(Q) while a^t^^ 
is 0(Q 5 ); thus, on the basis of scaling, cff£-body) ls expected to be much smaller than 
a \?>-body)- Meanwhile, the boost correction of Eq. (||) involves the expectation value inside 
the deuteron of two powers of momentum and a three-dimensional delta-function; and so 
scales as 0{Q 3 ). 

In Table ^ we exhibit various contributions for a characteristic wavefunction and their 
scalings in both baryon xPT and in the modified power-counting. The modified counting 
does a much better job of predicting the sizes of various mechanisms, presumably because 
it includes the full hierarchy of scales 7 -C M n <C A x , where A x sets the breakdown scale 
of the chiral expansion. 
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2a 
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0{Q A ) 


0.00184 


0.00175 


2bc 
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0(p 3 ) 


0(Q 5 ) 


-0.00017 


-0.00035 


^TodyV Ref - 35 1 


0{p 5 ) 


0(Q 5 ) 
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0.00041 



Table 4: Scaling of various contributions to the scattering length in baryon ;\PT and in 
the modified power-counting as discussed in the text. All matrix elements are evaluated 
with both the NLO (500 MeV) and the NLO (600 MeV) wavefunctions. 

In the modified power-counting we do not yet have the full result at 0(Q A ). In partic- 
ular, graphs of the type shown in Fig. f| are promoted from their 0(p 5 ) scaling in baryon 
%PT to 0{Q 3 ) in the modified counting. These graphs are corrections to a 2 ^_ body ^ with an 

insertion from at each vertex (Fig. f|(a)), or an insertion from C^L at one vertex and 
C^x at the other vertex (Fig. |(b,c)). Computation of these graphs is straightforward. 
Fig. |](a) evaluates to 

a %-W= 16 ^ +li /2) (^irt) (?) wf ' (19) 
while Fig. 0(b) and Fig. 0(c) are given by 



a (3-body) 



87^(1 + /X/2) 



4 5 + + 



32m 2 I 16vr 2 F2 




(20) 



The complete 0(Q 3 ) expression for the real part of the 7r-d scattering length is 

Rea nd = 2 J — ^2) fl+ a ( boost ) + a (3-6ody) + a (3-body) + a (3-body) (21) 
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where the various terms can be found in Eqs. (0), (|), (|12|) , ( |19"D and (|20D , respectively. 
The astute reader will have noticed that to the order we are working we can rewrite 
Eq. (^l|) in the form 



Re CL-n-a 



+ (1 + /i)f(a + ) 2 - 2(a-f 



2tt 2 



172 ) + a (boost) + 0(Q 4 



wf ■ 



(22) 



Thus, to 0(Q 3 ) the EFT with modified power-counting reproduces the simple intuitive 



picture of a single-scattering term plus a pion-exchange correction [|Tl[ . 

Notice that in the amplitudes of Fig. f|, the momentum dependence of the vertices 

(2) (3) 

from jQjy and has been neglected. While these corrections occur at the same order 
in baryon %PT as Eq. ([19]) and Eq. ( fZUD (O(jo 5 )), in the modified power-counting the 
momentum-dependent corrections begin at 0(Q 5 ) and are thus suppressed by two orders. 
These "range" corrections were recently computed in baryon xPT [|35|| . The size of these 
contributions is consistent with the 0(Q 5 ) scaling (see Table There are, of course, 
many other contributions at 0(Q 5 ), including two-pion-four-nucleon local operators with 
unknown coefficients Q. As we will see below, there are many more important contributions 
at 0(Q 4 ) which are easily computed in the EFT. Hence we will not consider any 0(Q 5 ) 
effects in isolation, but will instead perform the complete calculation up to 0(Q 4 ) in the 
modified power-counting. 



Figure 5: One-loop three-body diagram that contributes at 
0(Q 4 ) in the modified power-counting and at 0(p 5 ) in baryon 
XPT. 

It is surprisingly easy to compute the 0(Q 4 ) corrections to the 0(Q 3 ) result, Eq. (^2|). 
There are no further two-body effects from the 7r-N amplitude at 0(Q 4 ). There are, of 
course, diagrams similar to those in Fig. |4] and Fig. H(a), but with higher-order momentum- 
independent vertices. However, these corrections are already included in Eq. ( p2"D since 
they are subsumed in a ± . There is then a single three-body one-loop graph, Fig. |5|, 
representing the next term in the ttNN multiple-scattering series. This contributes at 
0(Q 4 ). Fig. |5| is one of many 0(p 5 ) baryon xPT graphs that contribute to a n d- However, 

4 It is, in a sense, unfortunate that these effects are suppressed, since one of the local operators at 0(Q 5 ) 
encodes information about the leading quark-mass dependence of the deuteron binding energy |16, 36). 
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the other graphs all have intermediate pions that are far from their mass-shell, and so are 
of higher order in the modified power-counting. 

We have evaluated Fig. |5] using the leading-order ;\;PT 7r-N vertex for ir-N scattering 
and have neglected contributions to the loop integral where pions are ~ M„ off their 
mass-shell. These neglected contributions are again higher order in the modified power- 
counting. This gives: 

a %-bodv) = C aam , ../o\ ( W7y?. I ( T^T ) • ( 23 ) 



64tt 4 (1 + (i/2) \2F% 



Notice that, as in Eq. (0), this matrix element has a Coulomb-like propagator, which 
signals the presence of an infrared logarithm. It is this enhancement that is captured by 
the modified power-counting. 

Our final formula for the 7r-d scattering length, valid to 0(Q A ) in the modified power- 
counting, is 



Rea - = 2 (TT^y ( a+ + (1+/i) 



[a + ? - 2(a-) 2 





+ a (boost) + 0(Q 5 ) . (24) 

where we have again subsumed higher-order effects into the 7r-N scattering lengths. Notice 
that since local two-pion-four-nucleon operators are not enhanced in the modified power- 
counting they appear at the same order as in baryon %PT, namely fifth order, and do not 
affect the fourth-order result given in Eq. ([24]) . 

Equation ([24]) is easily inverted to give a constraint on a + and a~ in terms of the 
experimentally-determined a n d from Eq. (f|). In order to evaluate the constraint we need 
the matrix elements (l/q 2 ) w f and (l/\q\) w f, as well as a^ boost ) (see Tables ^] and The 
boost effect does depend on the chiral low-energy constant, C2, but this counterterm was 
fixed in Ref. by fitting 7r-N phase shift data. The resulting relationship between a~ and 
a + is displayed in Fig. |5[. Notice that we choose to plot a + vs — a~ , in order to facilitate 
comparison with other extractions. 

The two bands in Fig. ^| correspond to two different evaluations of the deuteron matrix 
elements involved: one with the NLO wavefunction with an ultraviolet cutoff of 500 MeV 
(upper curve) and the other with a cutoff of 600 MeV (lower curve). In each case the finite 
width of the band arises from the error associated with the pionic-deuterium measurement 
of a n d. The spread in the values of C2 arising from Fits 1-3 of Table [L| also falls within 
these two bands. The result for a + found using "realistic" NN wavefunctions falls between 
these two bands, as, in most cases, would the effect of a typical 0(Q 5 ) correction. The 
experimental curves for the pionic- hydrogen shift and width are taken from Ref 0, while 
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Figure 6: Plot of a + vs —a~. The light shaded region and the 
dark band are from the experimental pionic-hydrogen width 
and shift, respectively, taken from Ref. ||. The dotted line 
encompasses the constraints from ir-N phase shift data and is 
taken from Ref. [[|. The dot is leading order %PT (current 
algebra). The two parallel bands are from Eq. (|4]) evalu- 
ated with the NLO wavef unction with an ultraviolet cutoff of 
500 MeV (upper curve) and 600 MeV (lower curve). 



the bounds from baryon xPT fits to 7r-N phase shift data are taken from Ref 0. From 
the overlap region in this plot we extract the S-wave 7r-N scattering lengths. We find 

a' = [0.0918(±0.0013) }M~ l ; 

a+ = [-0.0034(±0.0007) JM" 1 . (25) 

The theoretical error bars should be considered somewhat optimistic given our neglect of 
isospin violation. 



4 Comparison with other extractions 

The results of Eq. (^) are consistent with other recent extractions that do not make use 
of the EFT framework [[37], In particular, they are in agreement with very recent 
work of Ericson and collaborators |38| who find a~ = [0.0900(±0.0016) }M~ 1 and a + = 
[— 0.0017(±0.0010) }M~ 1 . The calculation of Ref. |38f includes the basic mechanisms 
present in the EFT employed here, but also includes additional contributions to which 
are higher order in the modified power-counting. 
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For instance, the analysis of Ref. [38 



includes a "dispersive contribution" (with a 
sizeable error bar) that is numerically comparable to our 0(Q 3 ) boost correction (but 
with opposite sign). This contribution accounts for most of the numerical difference 
between the two analyses. However, we find that the "dispersive contribution" occurs 
only beyond the order to which we have worked in the modified power- counting. 



I I 



Figure 7: Characteristic one-loop three-body diagram that 
contributes at 0(p 5 ) in x?T. This graph gives a contribution 
to the imaginary part of the n-d scattering length. 

One graph which contributes to both the dispersive and absorptive parts of a nc i is 
shown in Fig. [7|. This three-body one-loop graph is one of the leading contributions to 
Im a n d- It is 0(p 5 ) in baryon and 0(Q 6 ) in the modified power counting []. The 

0(Q 6 ) scaling of Fig. [7] is seemingly at odds with the experimental value of Im a n d given 
in Eq. (Q), which is consistent with an 0(Q 3 ) effect. It is easy to see why the scaling 
might break down for the imaginary part of the scattering length. The absorptive part of 
Fig. [7| is related by unitarity to the n d— > NN process at threshold: 



Im a„ 



— lim q a(ir 

47T I-*® 



d 



nn) 



(26) 



where q is the pion three-momentum in the 7r-d center-of-mass system. A study of tt 
d— > NN has been made in Ref. |39[| . Unfortunately baryon ;^PT applied to pion produc- 
tion significantly underpredicts the data. Presumably this is due to the large expansion 
parameter inherent in the pion-production process: since in this diagram momenta of 
order y/mM n are present, the baryon ;^PT expansion parameter is ^fp rather than p. 
Clearly this problem is present in any attempt to compute the absorptive contribution to 
the 7r-d scattering length in EFT, and is beyond the scope of this paper. Whether this is 
also an issue for the dispersive part is worth investigating, given the sizeable contribution 
claimed in Ref. 3~8|| . We should also mention that Refs. |37] and [[38] give an estimate of 
isospin violation that is consistent with an 0(Q 5 ) effect in the modified power counting. 
A systematic investigation of this very important issue within EFT is needed. Finally, 
our result for the isoscalar scattering length is consistent with a previous baryon %PT 
determination at 0(p 3 ) | 21| . Because a hybrid approach employing phenomenological 
wavefunctions was used in that paper, no theoretical error bar was given. 

5 Since pion number is not conserved at the interaction vertices, this graph appears as a two-pion-four- 
nucleon local operator in H7rEFT. 
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5 Conclusion 



In this paper we have developed an EFT for the 7r-d scattering length. While naively, 
baryon ;\PT, as modified by Weinberg jn] for the NN system, captures the correct hi- 
erarchy of scales, we have seen that there are several subtleties in the analysis of 7r-d 
scattering near threshold, which suggest a more practical power-counting. The new in- 
gredient is the observation that one may develop a hierarchy between the pion mass, 
M n and the deuteron binding momentum, 7. One can then pick out the pieces of the 
amplitude that become large in the formal limit in which 7 goes to zero and M n is held 
fixed P^fl . This does not, of course, indicate that anything is amiss with baryon %PT 
power-counting. Baryon %PT should work just fine for n-d scattering in the threshold 
region, but high-order calculations are required if truly accurate results are desired. The 
advantage of the modified power-counting, and of EFT generally, is that it immediately 
isolates the large contributions, without requiring explicit calculations of matrix elements 
that end up being smaller than the theoretical error. 
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